Physics on the edge: contour dynamics, waves and solitons in the quantum Hall effect 
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We present a theoretical study of the excitations on the edge of a two-dimensional electron 
system in a perpendicular magnetic field in terms of a contour dynamics formalism. In particular, 
we focus on edge excitations in the quantum Hall effect. Beyond the usual linear approximation, a 
non-linear analysis of the shape deformations of an incompressible droplet yields soliton solutions 
which correspond to shapes that propagate without distortion. A perturbative analysis is used and 
the results are compared to analogous systems, like vortex patches in ideal hydrodynamics. Under a 
local induction approximation we find that the contour dynamics is described by a non-linear partial 
differential equation for the curvature: the modified Korteweg-de Vries equation. 
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I. INTRODUCTION 

The theoretical description of many-body systems is 
often best realized in terms of collective modes, i.e. the 
familiar sound waves in solids or plasmons in charged 
systems. Collective modes become especially important 
when their energies are lower than competing single- 
particle excitations. Sometimes, however, both single- 
particle and collective modes in the bulk of a system are 
gapped or scarce and these systems are often referred to 
as "incompressible." This incompressibility can be real 
or a convenient limit due to large differences in relevant 
length- or time-scales, as in the case of the macroscopic 
motion of a liquid. Under these conditions, one can usu- 
ally focus our attention on the motion of the boundaries 
of the system, which will generally have softer modes, 
with frequencies much lower than those in the bulk (e.g. 
surface waves in a liquid droplet travel at speeds consid- 
erably slower than sound waves). 

Concentrating on the motion of the boundary of the 
system has a considerable advantage: the reduction in 
the dimensionality of the problem often permits simpler 
analytical treatment, or a tremendous reduction in the 
effort needed to numerically solve or simulate the prob- 
lem. Associated with this incompressibility, however, one 
usually finds microscopic conservation laws that translate 
into global constraints on the whole system, even when 
the microscopic dynamics is completely local (e.g. the 
volume of a liquid droplet is conserved). These global 
constraints enter the edge dynamics through Lagrange 
multipliers associated with the conserved quantities .EJ 
These conservation laws are often evident when the mo- 
tion of the system is observed, and it is interesting to see 
how they are embedded in the dynamics of the bound- 
ary; that is, how these essential aspects of the problem 
are related to the laws of motion of the edges. 

These shape deformations, and their dynamics, have 
played an important role in the understanding of numer- 
ous problems in diverse fields of physics. The incompress- 
ibility is reflected in the existence of a field which is piece- 
wise constant, so that there is a sharp boundary between 



two or more distinct regions of space with different physi- 
cal properties. This field can be of classical origin like the 
density of a liquid or the charge density of a plasma, or it 
can originate in the quantum mechanical properties of the 
system, like the magnetization of a type-II superconduc- 
tor. There are various examples where these questic 
are relevant, such as waves on the surface of a liquid, E 
the motion of non- neutral plasmas Jj low lying "rotation- 
vibration" modes of deformed nuclei,cl the evolution of 
atmospheric plasma clouds ,13 pattern formation an— ferro- 
magnetic fluids Jj vortex patches in ideal fluids JjEj and 
two-dimensional electron systems (2DES) in strong mag- 
netic fields £3x3 

The edges of a two-dimensional electron system 
(2DES), and in particular the edges of a quantum Hall 
(QH) liquid, present a unique opportunity to study the 
dynamics of shape deformations in a clean and controlled 
environment. The 2DES in the QH state is incompress- 
ible, so that the electron density is approximately piece- 
wise constant, suggesting that a contour dynamics ap- 
proach to studying the droplet excitations is viable. The 
lack of low-lying excitations eliminates dissipative effects, 
further simplifying the treatment of the problem. In ad- 
dition, the charged nature of the system facilitates the 
excitation and detection of deformations of the droplet. 

In this paper we will formulate the study of the excita- 
tions of a droplet in a 2DES as a problem in contour chu. 
namics. In the usual treatment of the edge excitations ,Ej 
a linearization of the equation of motion is done at early 
stages, thus limiting the applicability to small deforma- 
tions of the edge of the system from an unperturbed 
state. In this paper we consider non-linear terms which 
are present in the full contour dynamics treatment. We 
first present perturbative results for non-linear deforma- 
tions of the 2DES shape. For the sake of comparison, 
and as means of verification, we also apply this method 
to the vortex patch case, which has well known exactEl 
and numericallill solutions. We then show that the lo- 
cal induction approximation to the full contour dynam= 
ics generates the modified Korteweg-de Vries (mKdV)Ej 
equation for the curvature dynamics; the mKdV equa- 
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tion also arises in studies of vortex patched and sus- 
pended liquid droplets. □ The mKdV dynamics conserve 
an infinite number of quantities, including the area^-aen- 
tcr of mass, and angular momentum of the dropletjlj so 
that our local approximation to the nonlocal dynamics 
preserves the important conservation laws. The mKdV 
equation also possesses soliton solutions, including trav- 
eling wave solutions. 

In Sec. [n] we present a brief review of the hydrody- 
namics of a two-dimensional electron system in a strong 
perpendicular magnetic field and analyze the bulk and 
edge excitations in the linear approximation. Sections 
III and IV and analyze in more detail the dynamics and 



kinematics of these edge modes, and ask what conditions 
must be satisfied so that a large edge deformation is able 
to travel without any dispersion, that is, preserving its 
shape. The question is posed in terms of a non-linear 
eigenvalue problem and is solvedperturbatively to fifth 
order in the deformation in Sec. |\] (for completeness, we 
also solve the analogous problem of vortex patch defor- 
mations in Appendix |b]). Some solutions and limiting 



cases are then presented in Sec. VI. An alternative ap- 



proach to find non-dispersive or invariant deformations 
of the edge, namel y th e local induction approximation, is 
developed in Sec. VII , where results are also compared 



with the perturbative approach of Sees. and VI 
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that is, bulk magnetoplasmons are gapped, and in the QH 
regime can be neglected since hui c ^ 17meV~200K at B 
= 10 T (in GaAs), whereas T ~ IK. We will disregard 
them from now on and concentrate on the excitations at 
the edge which, as we shall see, have considerably lower 
energies. 

The theory of small deformations of the edge has been 
extensively studied.EZrE] The main conclusion is that for 
strong magnetic fields, when Landau level quantization 
becomes important, the only low- lying modes are edge 
modes which propagate in only one direction along the 
edge of the 2DES. We further classify these modes into 
the "conventional" edge magpetoplasmon mode, with a 
singular dispersion relation:H3 
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where k is the mode wave-number, 7 sa 0.5772 is the Eulcr 
constant, and a is a short-distance cut-off (the largest of 
the transverse width of the 2DES, the magnetic length, 
or the width of the compressible edge-channel). In addi- 
tion, for wide compressible edges, "acoustic" modes can 
be found: 



II. HYDRODYNAMICS OF A 
TWO-DIMENSIONAL ELECTRON SYSTEM 

Consider a 2DES in a perpendicular magnetic field. 
Treated as a classical fluid, the system is characterized 
by the electron density n(r) and velocity v(r). If we 
neglect dissipative effects (which is essentially correct in 
the quantum Hall regime), the dynamics is determined 
by the Euler and continuity equations: 
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where oj c — eB/m e is the cyclotron frequency and e is 
the dielectric constant of the medium. The first term 
on the right-hand side of Eq. (|l|) represents the Lorentz 
force, the second is the Coulomb interaction, and E ext is 
the electric field due to the background positive charge, 
gates, etc. 

Consider first the possible bulk collective excitations 
in a uniform system. These are oscillations of the den- 
sity around the equilibrium solution v = 0, n — h. If we 
consider small amplitudes and Fourier analyze Eqs. (|l|) 
and (0), we find that 



ik-r— iujt 



v(r,t) = voe 

Ti 

<5n(r, t) = — k-vo e l 

UJ 



(3) 



2 fie 2 
/ 

em e uj c j 



J 



1,2,. 



(6) 



These results are approximately correct as long as iner- 
tial and confining terms are negligible. The first requires 
that ne 2 /em e a sSj w^, which is usually true in quantum 
Hall conditions .Ej In addition, the effects of the confining 
potential have been neglected. While the confining po- 
tential is essential for long-term stability and is usually 
not negligible when compared to interactions, its effect 
is mainly reduced, in the simplest cases, to an additional 
i> cx t k term in the frequencies, where the external velocity 
is given by 



m e uj c 



■e,xE, 
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Recent time-of-fiight measurements^!] in 2DES confirm 
this simple picture. The theoretical formulation above, 
however, is restricted to small deformations of the bound- 
ary. In what follows we will consider a formulation that 
goes beyond this limitation. 



III. DYNAMICS OF THE EDGE MODES 

Consider the case when the edge between bulk and 
vacuum is sharp. The electronic density has some con- 
stant value n inside the edge and vanishes outside (see 
Fig. [I]). Since the edge is essentially one-dimensional 
only the edge magnetoplasmon is important (and even 
in more general cases this mode is the most readily 
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observablal3'EJ). -Following the derivation of the edge 
magnetoplasmonai3c3 we neglect inertial terms in Eq. 
(0), thus obtaining an equation for the electron velocity: 



v(r) = - 



V x e z / d z r 

em e LO c J A |r - r 



2j n ( r ') 
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where A is the area of the droplet (see Fig. [I]). 



equilibrium 
deformed shape 




n=0 



FIG. 1. A charged incompressible liquid in a magnetic 
field. We assume a piecewise constant electron density {n = n 
inside, while n — outside the droplet). The parameteriza- 
tion R(ip,t), the tangent t and normal n unit vectors to the 
boundary F are indicated; s is the arc-length and 8 the tan- 
gent angle. 

Let us now concentrate on the "internal" velocity given 
by the first term in Eq. (||). We should note that the ne- 
glected external velocity v ex( derived from the confining 
field [Eq. (^|)] is important for long term stability of the 
system, and will modify the propagation velocity. In gen- 
eral, it could also change the shape of the modes, yet one 
can devise situations in which this effect is not important, 
i.e. a linear confining potential in a rectilinear infinite 
edge, or a parabolic potential for a circular droplet. 

For an incompressible 2DES with a piecewise constant 
density, the density can be taken outside the integral; 
then using Stokes' theorem, the area integral can be 
transformed into a line integral over the boundary T = OA 
of the electron liquid: 



v(r) 



ds 



, too 



em e u c J T |r - r(s')| ' 
Here t(s') is the unit tangent vector at arc-length s': 



(9) 



t(s') = d s ,r(s') = iv 
n(s') = -e z x t(s') , 



(10) 



and we defined the unit normal vector n(s') for later use. 
The short distance singularity in the integrand is cut off 
at a length scale r Q . Equation (^|) forms the basis of our 
contour dynamics treatment — it expresses the velocity of 
the edge in terms of a nonlocal self-interaction of the 
edge. 

We see that: (i) the dynamics is chiral, being deter- 
mined by the tangent vector; (ii) the fluid contained 



within r is incompressible, so that the area is conserved. 
It is interesting to note similar traits wjerp-found in the 
past for the dynamics of vortex patches;Eaii3 indeed, it is 
this analogy which inspired the present work. A detailed 
description of the vortex patch case, and stressing the 
similarities and differences with the shape deformations 
of the 2DES is presented for completeness in Appendix 



IV. KINEMATICS OF UNIFORMLY 
PROPAGATING DEFORMATIONS 

Having determined the velocity of the electron liquid, 
we now focus on the motion of the 2DES boundary T. 
The velocity of a point on the boundary can be written 
in terms of the normal and tangential components, 



U{s)n + W(s)t. 



(11) 



The tangential velocity W(s) is largely irrelevant, as it 
solely accounts for a reparametrization of the curve; the 
boundary motion is determined by the normal compo- 
nent of the velocity. 




: 
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FIG. 2. A uniformly propagating edge deformation. The 
boundary satisfies R(<p,t) = R(ip — fit). 

We now ask whether there are modes which propagate 
along—the boundary with no change in shape. Previous 
workEO has focused on small perturbations of a straight, 
infinite edge. Here we consider deformations of a circular 
droplet of incompressible electrons (the non-linear defor- 
mation of a straight infinite edge is discussed in Appendix 
[f]) . A uniformly propagating mode is, in this case, char- 
acterized by a boundary that moves like a rotating rigid 
body (see Fig. ||), namely, the radius of the boundary 
satisfies 



R(<p,t) =R(<p-ttt), 



(12) 



where ip is the azimuthal angle, and f2 is the angular fre- 
quency of the boundary rotation. This translates into a 
condition for the normal velocity: 

U = n(r) • v(r)| rer = fi n(r) ■ {e z x r) . (13) 

We seek surface shapes R(ip) (see Figs. |l| and |^) that 
rotate uniformly, satisfying Eq. (|l3|). Consider the fol- 
lowing parameterization of the surface: 



R(<p) = Ro[l 
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In this parametrization, we can write the unit tangent 
vector explicitly as 



t(¥>) = 



where 



T { V ) = ^^- =er#(<p)+e v R(<p). 
Likewise, the unit normal vector is given by 



n(tp) = e z x 



■(v)l 



(15) 



(16) 



(17) 



Given the identity r s ds = r v dip, the normal velocity of 
the boundary, derived from Eq. (0), can be written as 



em e uj c J |R(<p) -R(</)| 

while inserting Eqs. © and @ into Eq. © yields 

(19) 



{/(<?) = n(R) • (e 2 x R) 



-Ml 



I 
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We seek frequencies ft and coefficients bi that satisfy Eqs. 
( |l8| ) and ([l9]). The solutions to this non-linear eigenvalue 
problem represent edge modes that propagate without 
any dispersion. 



V. SOLVING THE NON-LINEAR EIGENVALUE PROBLEM— PERTURBATIVE APPROACH 

Unfortunately, it has not been possible to solve exactly the non-linear eigenvalue problem for bi and 91 as written in 
Eqs. ( |l8| ) and (|l9|). We therefore seek a perturbative solution by expanding the right-hand side of Eq. ( |l8|) in powers 
of b^ This allows us to to go beyond the linear approximations used in the past, and we have succeeded in calculating 
shape deformations to 0[bf] and angular frequencies to C[6f]. Expanding the non-linear eigenvalue problem to fifth 
order, we find the condition 

& f h + - ^2 °i-pbp J = Ql h + ^2 i?;-p, p h-pbp +^ S , j_ P)P _g )9 bi^ p b p ^ q b q 
V p / p p,q 

~l~ ^ T-'l—p^p—q^q- J r.rbl~pbp^qbq— r b r -\- ^ ^ p,p— q.q-^r.r— s,s^/— pbp~ qbq— rb r — s b s -\- (D[b^ , (20) 



p,q,r 



p,q,r,s 



where 



(21) 



and the "matrix elements" Q, i?, 5, T and [/ are obtained from Eq. ( |lS| ) in Appendix [A|. 

Since the equation to be solved results from an expansion in powers of it is sufficient to solve Eq. ( pp| ) pertur- 
batively, by expanding in powers of the largest coefficient bh- Let us assume that 



= 0[8], for l = ±L. 
= 0[S 2 }, for l^±L. 



We then consider an expansion of the coefficients bi and eigenvalue 91 in powers of 5: 
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b, = I \ ' 



91, 



n 



(2) , 7,(3) 
(0) 
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for l = ±L. 
for l^±L. 
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(i) 



(2) 



(22) 

(23) 
(24) 



where 6( fe ) and 91^ are of order 0[<5 fe ]. By substituting expressions ( p3| ) and ( p4[ ) into Eq. j20|), and grouping terms 
order by order, we find that: 



1. The first-order term yields the dispersion relation in the linear approximation: 

9l^[L] = Q L : 



(25) 



4 



2. Second-order terms couple the fundamental and second harmonics of the deformation, with no correction to the 
eigenvalue: 



QV>[L] = 0, 



€\ L = i 2 * L : L ^ 



Ql 

2 Ql — Q2L 
b\ 2) = , for I ^ ±2L : 



(26) 



3. Third-order terms give the first correction to the frequency and couple the first and third harmonics: 

f2 (2) [L] = ZSl^-l^l? + (2R 2L ,-l - QL)bf L , 



41 - 

6f) = 0. 



Ql — Qzl 
for I ^ ±3L ; 



(27) 



4. Fourth-order terms provide coupling to both second and fourth harmonics: 
ft (3) [L] =0, 



b±l L = 1 

Ql — Q2L 



(2i?; 



3L,-L 



Q L )bf L b^-^[L] + 



6 (4) _ 



^ {i(2i?2L,2L - Ql)^] 2 + (2i?3LX - QzO&S&j? + 3S 3i ,z, L &g [b^] 2 + T L , L , L , L [fc^] 4 



(28) 



6, (4) = , for I ± ±4L ; 



5. The fifth-order term follows the same pattern. The couplings to higher harmonics is quite complex and we chose 
not to pursue it. We show only the correction to the eigenvalue: 



n (4) [I] = (2R 2 L,-L - Ql)4l + {2RiL,-2L - Ql 



6 (3)j,(2) 
U 3L u 2L 

1.(1) 



+65 l ,2L,-2l[6 2 l] 2 



(2) 
2L 

■ 2 ) r^C^is 



■5*31, -L ,-l0 3L o l 



^(1)14 



(29) 



Note that the largest term in the perturbative expansion corresponds to the lowest harmonic and that higher order 
elements preserve the rotational symmetry Cl (rotations by 2n/L) of the initial term. 



VI. SOLUTIONS FOR THE TWO-DIMENSIONAL 
ELECTRON SYSTEM 

We now show some invariant shapes for the 2DES (see 
Sees. 



Ill 



and IV). Coefficients 6/ and frequencies f2 are 
obtained using the formulae of Sec. and matrix elements 
Q, R, S, T and U calculated in Appendix [A|. 

Table | summarizes these results for states of rotational 
symmetry Cl, with L = 2, 3, 4 and 5. Some particular 
shapes obtained from these coefficients and Eq. (14) are 
shown in Fig. I [for a comparison with similar states for 
vortex patches see Figs. I and § . 

For large deformations, the appearance of oscillations 
indicate that higher order terms are needed, since the 
perturbative method corresponds to a truncation of the 
Fourier decomposition [Eq. (p"4|)]. An alternative ap- 
proach, based on a local induction approximation, pro- 



vides a better description in those cases. This alternative 



formulation is presented in Sec. VII 



It is interesting to note that the linear result for the 
frequency is [see Eqs. (p5|) and (A3)] 



nW[i] = Q l = 4j2 



1 



k=2 



2k- 1 



(30) 



where ths last sum can be related to the digamma 
functioned ijj{\l\ + l/2). This lineax.result has been previ- 
ously derived by several authors;E3 corrections are 0[bf}. 
For a direct comparison with Eq. (^J), which corresponds 
to the large-l limit, we substitute the asymptotic expan- 
sion for the sum in Eq. (30); multiplication by Rq yields 
the propagation velocity: 



5 



du>o{k)/dk obtained from Eq. (|5|) after the substitu- 
tion I ~ ka. The dispersion for these linear edge exci- 
tations has, been confirmed, experimentally in both the 
frequencjcS and timeEd'oEa domains. 



TABLE I. Stationary deformations of a two-dimensional electron system. Angular frequencies and lowest harmonics of the 
deformation for L = 2, 3, 4, 5. 
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Q. = (em e u) c Ro/ne 2 ) ft 


b%i 


63; 


64; 


2 
3 
4 
5 


| - 5.09048 b\ - 16.1333 6| 
f ~ 16.8464 b% + 39.6332 bj 
f§| - 30.0989 b\ - 154.821 b\ 
§f| - 45.5507 b\ - 418.524 b% 


§§ 6| + 1.86911 b\ 
1.63473 b% +8.36452 b\ 
1.94074 b 2 4 + 24.3045 6| 
2.17945 6| + 50. 1748 6| 


0.911769 6| 
1.97442 6i 
2.56102 6| 
2.74603 6| 


-0.345388 b\ 
-0.540557 b\ 
-3.45272 6t 
-9.80691 6| 
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which closely corresponds to the group velocity w„ 



VII. LOCAL INDUCTION APPROXIMATION 

As we have seen, the motion of the edge is determined 
by the velocity of the fluid at the surface. The nonlocal 
equation for the velocity of the boundary, Eq. (^), can 
be turned into a differential equation for the curvature of 
the boundary if we concentrate on the locaLcontributions. 
This local induction approximation (LIA)c3 was explored 
by Goldstein and Petrich in a scrics-of-papers dealing 
with the evolution of vortex patches S3t3 The situation 
is considerably more favorable in this problem, due to 
the more rapid decay of the interaction [1/r for charges 
vs. ln(r) for vortices, see Eqs. @ and (B3)]. Figure [j] 
defines most terms used in this section. 

The LIA corresponds to the introduction of a large dis- 
tance cut-off A in the expression for the velocity of the 
boundary v[r(s)] [Eq. ([}])]: 



{■■.}ds' 



s+A/2 

{■ 

-A/2 



•} ds' . 



(32) 



The line integral in Eq. (||) is then calculated by ex- 
panding the integrand in powers of A = (s'— s). By using 
the Frcnct-Scrrct relations 



= t, 
= — k n , 
= nt , 



(33) 



where k = 1 
obtain 



is the local curvature of the boundary, we 



t( S ')^t( S ) 

+n(s) 



A 2 2 A 3 



(34) 



A 2 

■A/! K, 

2 



r(s') - r(s) ~ t(s) 



-n(s) 



A A 3 2 A 4 
A 3 



(35) 



A 2 _ _ ; 



To lowest order the normal and tangential velocities 
then are given by: 



U UA = 

W L ia = 



em e Lo c 



A 2 



2r 



11A 2 
96 



(36) 



It is worth noting that since the rate of cha nge of the area 
A of the droplet \sA t = § ds U{s) [see Eq. (jof)], the LIA 
with Uua oc k s (or any exact differential) automatically 
conserves area. This is not surprising since we started 
with an incompressible system, but shows that the local 
approximation used has not introduced an obvious error. 
It is also interesting to realize that the perimeter C of the 
curve derived from these velocities is conserved as well: 
C t =§ ds (kU + W s ) = [see Eq. (|D7|)]. 

The time evolution of a curve in two dimensions 
is given— .ouite generally by the integro-differential 
equationlla-EZl (see Appendix [5]) 



ft'/ 



[k 2 + d ss ]U + k s W - k s [ [kU + W s ,]ds' . 
Jo 



(37) 
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We now introduce the results from Eq. (pq). A "gauge" 
change, realized by modifying W so that the integrand 
of Eq. ( [37] ) vanishes, eliminates the remaining non-local 
dependences and yields 
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8 \ 



3 2 

.2 s 



(38) 



After a simple ti, 
mKdV equation 



rescaling, the curvature satisfies the 



(39) 



The mKdV dynamics are integrable, with an infinite 
number of globally conserved geometric quantities,t3 the 
most important of which are the center of mass, area, 
and angular momentum of the droplet. 

The mKdV equation possesses a variety of soliton solu- 
tions, including traveling wave solutions and propagating 
"breather" solitons (see Appendix |]) . Here we will focus 
on the traveling wave solutions of Eq. ([39]) of the form 

k(s, t) — k(z) , with z = s — ct , (40) 

which represents uniformly rotating deformed droplets 
(see Fig. ||). The ordinary differential equation for k(z) 
can be integrated twice with the result 



^ 2 = -l« 4 



26, 



(41) 



where a and 6 are constants of integration (a = 6 = for 
infinite systems, see Appendix |]). This ordinary differ- 
ential equation can be easily integrated: 



V(k) 



= ± 
4 



C K 



d K ' 

2a n + 46 . 



(42) 
(43) 





t = 



t>0 



FIG. 3. A uniformly propagating edge deformation. The 
curvature satisfies K,(s,t) = k(z), with z = s — ct. 

This problem is thus analogous to a particle moving in 

Involved can 
_l and depend 



a quartic potential (Fig. |4]). The integra 
be expressed in terms of elliptic integral: 
crucially on the nature and location of the zeros of V(k). 
For curves with finite perimeter and no self-crossings we 
find that it is necessary to have two real and two com- 
plex zeros: K max , K min , and -(K max + K min )/2±i£ (see 
Appendix |e|). Note that K max and K m i n correspond to 
the maximum and minimum curvatures of the boundary. 




FIG. 4. Some possible potentials V(n) and corresponding 
schematic solutions for the shape of the boundary. Potentials 
with two real zeros may have physical solutions: (a) Curva- 
ture is fixed, the solution is a circle; (b) Curvature is positive 
definite, the boundary is convex; (c) Positive and negative 
curvatures, (c) Potentials with four zeros have solutions with 
unphysical self-intersections. 

The periodic solutions of-J3q. ([42]), expressed in terms 
of Jacobi elliptic functions ,E3 are given by (Appendix |j): 



k{z) 



where 



(Q^max+PKmin) + (P^min — 

gK max )cn( % /Mf I A) 



(p + q) + (p- q) cn( % /Mf | A) 



(44) 



p = 


\/(3/^max 


H~ ^min) 2 " 


-e, 


q = 




- 3K m i n ) 2 - 


-e, 


A = 






(p - q) 2 ] l^pq 



(45) 



The free parameter £ is actually determined by the 
boundary conditions. The period of k is given by the 
elliptic integral 



L K = 



pq 



K(X). 



(46) 



We now require that the tangent angle increases by a 
factor of 27r after an integer multiple / curvature periods, 
so that the curve is closed and with no self crossings (see 
Fig. |: 



0(lL K ) 



IL. 



k(s) ds = 2-7T . 



(47) 



It is evident that the resulting curves are invariant under 
rotations by 2n/l, that is the curves have Ci symme- 
try. The curves thus generated can be characterized by 
(I, K max , K m in) or more conveniently, although indirectly, 
by the symmetry, the area and the perimeter of the curve. 
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FIG. 5. Illustration of the boundary conditions implied 
by Eq. (|4T[) . Top: after I periods L K of the curvature k(s), 
the tangent angle 9 increases by precisely 2n. Bottom: the 
resulting deformed boundary. 

The contour shape can be easily determined once the 
tangent angle 9(s) is known as a function of arc-length. 
Since dz = dx + i dy = exp[z 0(s)] ds, we have 

z(s) = x(s)+iy(s)= f e ie(s 'Us' . (48) 
Jo 



The full lines in Fig. |6| show some uniformly rotating 
soliton shapes, calculated from Eqs. ( ft4| ) and (|48|) . These, 
are essentially identical with Goldstein and Petrich'sll3 
soliton solutions for the vortex patch problem, but the 
more local nature of the interaction in the 2DES case 
should guarantee a better correspondence with the exact 
solutions including non-local terms. Indeed, the curves 
resulting from the perturbative method and the LIA are 
quite close, even for considerable deformations of the 
boundary. For larger deformations the perturbative re- 
sults show artifacts due to the limited number of Fourier 
components. The advantage of the LIA becomes evident 
in this case, since it is an expansion in powers of the 
curvature and not the deformation, and thus allows for 
relatively large long- wavelength deformations. More sig- 
nificantly, the LIA and the resulting integrable dynam- 
ics allow one to uncover geometrical conservation laws, 
which would be hidden in a perturbative calculation.Ea 
This advantage comes at a price: the detailed informa- 
tion on frequencies is obscured by the introduction of the 
long distance cut-off A and by the gauge transformation 
of the tangential velocity W, while the frequency is easily 
obtained in the perturbative calculation. 




FIG. 6. Uniformly rotating shapes of a 2DES. Solid lines: solutions of the mKdV equation obtained from the local induction 
approximation. Dotted lines: solutions obtained using the perturbative expansion. The values of the coefficient and the 
ratio of curvatures cr = «: m i n /K max are: (a) b 2 =0.073, a = 0.4; (b) b 2 =0.19, a = 0; (c) 6 2 =0.36, cr = -0.2; (d) 63 = 0.027, a = 0.4; 
(e) 63=0.10, cr = -0.2; (f) fe 3 = 0.29, cr = -0.45; (g) b 4 = 0.014, a = 0.4; (h) fe 4 = 0.089, cr = -0.4; (i) &4 = 0.24, cr = -0.56. 
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VIII. CONCLUSIONS 

A contour dynamics formulation of the excitations on 
the edge of a two-dimensional electron system in a mag- 
netic field has allowed us to demonstrate the existence, 
beyond the usual linear regime, of shape deformations 
that propagate uniformly. A local approximation to the 
nonlocal dynamics shows that the curvature of the edge of 
the droplet obeys the modified Korteweg-de Vries equa- 
tion, which hasjjjitegrable dynamics and soliton solutions. 
Earlier studiesEa of edge channels in QH samples have 
shown the presence of nonlinear waves, but in Ref . ^9] the 
nonlinearity originates in the variations of the strength of 
the confining field (the E cx t of Sec. ||), whereas here we 
concentrate on nonlinear effects originating in geometri- 
cal effects. 

Since these solutions are dispersionless, it may be pos- 
sible to distinguish them from linear edge waves in timc- 
of-flight measurements. of the type depicted in Fig. [?]. In 
a circular QH systeiro a voltage pulse applied to a gate 
produces an edge deformation. The deformation propa- 
gates along the edge of the system in one direction and 
is detected, i.e. by means of a capacitive probe. For 
this geometry, the pulse comes back repeatedly and it 
is possible to observe it after numerous passes. While 
a gradual decrease of the amplitude is always expected 
due to residual dissipation, dispersive and non-dispersive 
modes may be distinguished by the preservation of the 
general shape and width of the non-dispersive modes. 

On the theoretical side, it would be interesting to con- 



nect our hydrodynamic treatment of these edge solitons, 
with field-theoretical treatments of edge excitations Ii3££l 




Vout dispersive 




FIG. 7. Soliton excitations travel without shape changes. 
In a circular QH systemcJ with non-contact probes it is pos- 
sible to observe repeated signals. While a gradual decrease 
of the amplitude is always expected due to residual dissi- 
pation, dispersive and non-dispersive modes may be distin- 
guished by the preservation of the general shape and width of 
the non-dispersive modes. 
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APPENDIX A: MATRIX ELEMENTS FOR THE NON-LINEAR EIGENVALUE PROBLEM 



In order to determine the matrix elements Q, R, S, T and U for the non- linear eigenvalue problem [Eq. (p0[)], we 
need to expand Eq. (|l8| ) in powers of the coefficients bi of the parametrization [Eq. (|lj)] . We initially realize that 



n(ip)-T(<p + uj) = 



— sin lu 



2^6/e 



iltp 



Ilu IlU 
cos — sin lu + I cos uj sin — 
2 2 



2 x (l-2p)uj . „ . . (l-2p)w 
(lp—p —1) cos smu)~{l — 2p) coslu sin 



(Al) 



|R(^ + c)-R(^ 



4 sin^ — 



1 + 2 ]T 6, e'^ 2 cosf-£ b^b p e^/ 2 — 1 ^— JL- 



. (l-p-l)u ■ (p+l) 

sin i — F 2 1 sin KF ' 



2 sin 2 f 



(A2) 



By expanding xi{(p)-T(tp + o;)/|R((p + lu) — H((p)\ in powers of bi, and integrating over lu [see Eq. (fL8|)l, one obtains 
the relevant matrix elements Q, R, S, T and U. All integrals converge without the need for short distance cut-offs. 
The three lowest order terms can be written as follows: 



(A3) 



k=2 







w 1 \i-p\ 1 \p\ 1 



(A4) 



S, 



!-p,P-9,9 



fe=l 

■5 

7 
1 

12 



fc=i 



fc=i 
'7 



Z-p 



i-3 



p-g 



l-p + q 



1-il 2 l-4p 2 1-V l-4(Z-p) 2 l-4(Z-g) 2 l-4(p-g) 2 l™4(Z-p + g ) 2 

-(^« 2 )i:^-< i+ * 2) i:^ + < 5+4 ' J >i:2^ + [ 5+4 <'-^ l 2 l 2^ 



fc=i 



fc=i 



ip-?i 



ij-p+«i 



[i+4(/- (Z ) 2 ]^^_ + [5 + 4(p- (7 ) 2 ]^^_-[i+4(/-p +9 ) 2 ] y: 



fe=i 



Higher order terms are long and uninspiring. 

APPENDIX B: COMPARISON WITH THE 
VORTEX PATCH CASE 

For the sake of comparison, we draw analogy to the 
case of a vortex patch, a two-dimensional, bounded re- 
gion of constant vorticity uj p surrounded by an irrota- 
tional fluid.E The vorticity can either be distributed, as 
in a regular fluid, or concentrated in individual vortices, 
in the case of a superfluid (in this case it is clear that the 
hydrodynamic treatment will be valid only for length- 
scales larger than the inter- vortex spacing). The impor- 
tant thing is that in ideal fluids the area of the vortex 
patch is conserved due to Kelvin's circulation theorcm,c3 
and therefore the vortex patch is essentially incompress- 
ible. Figure || can be used to describe this case by re- 
placing the electron density n by the vorticity ui p . 

The steady state solution in this case is clearly a cir- 
cle, and small deformations of the boundary travel along 
the boundary itself, as has been known for a long time.Ll 
One can also ask what happens when the deformations 
are large: are there modes that do not change shape, i.e., 
solitons? One such solution has been known since Kirch- 
hoff 's time: an ellipse with constant vorticity will rotate 
uniformly in an ideaJpfluid. Numerical calculations by 
Deem and ZabuskyE3llil in the 70's obtained additional 
invariant shapes. 

For inviscid incompressible fluids, the equations of mo- 
tion for the fluid are simply given by: 



k=l 



k=l 



V • v = 0, 
Vxv = w, 



(Bl) 
(B2) 



,(r) = -^ (fds't(s') In 



|r - r(s') 
ro 



(B3) 



(A5) 



where u> is the vorticity. Assuming that the velocity far 
from the patch vanishes, the velocity of the fluid in pres- 
ence of a region of finite vorticity lo p can be expressed 
as 



The arguments presented in Sees. [II and IV are now 
applied mutatis mutandis to this case. The only differ- 
ence with the 2DES comes from the fact that the kernel 



in the interaction is now logarithmic, and Eq. (Q) is then 
replaced by Eq. (B3). We see that, as in the 2DES case 
(Sees. Ill and IV), the dynamics is chiral, being deter- 
mined by the tangent vector; and since the fluid con- 
tained within r is incompressible the area is conserved. 
The normal velocity of the contour is then given by 



27T Jo 



|R(¥>)-RV)I 



(B4) 

As before, we seek solutions that satisfy Eqs. ([TJ|) [or 
@] and @. First, we follow the procedure of Ap- 
pendix [X] to determine the matrix elements Q, R, S, T 
and U. The first few of these are given by 



oV 
l-^p,p—q,q 



1 l 

2 ~ 2|7| 
1 

~W\' 

1 r 



(B5) 
(B6) 
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1 



\2l-q\-\p- 



-p\ + \2l-p\ 

q\ + \i + p- 



+6(2l-p) + S{l-q)+S{l+p-q) . (B7) 

We then apply the results of Sec. |V|, to determine the 
amplitude of the lowest harmonics for the vortex-patch 
case, using the appropriate matrix elements [Eqs. ( |b!j| ) 
(B7)]. These results are summarized in Table |n|, and 
some of the resulting invariant shapes are shown in Figs. 
| and §. „ 

It has long been knowro that an elliptical region of vor- 
ticity in an otherwise irrotational fluid, namely the Kirch- 
hoff ellipse, rotates uniformly with angular frequency 
fl = (wj,/4)(l — b 2 /a 2 ), where (a + b) and (a — 6) repre- 
sent the maximum and minimum radii respectively (see 
Appendix ^). A simple analysis shows that both fre- 
quency and angular components bi shown in Table || 
exactly match a series expansion of the ellipse. Figure 
H compares the perturbative results with the exact solu- 
tion. Even for relatively high deformations both results 
are in reasonable agreement. 
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FIG. 8. The I = 2 stationary deformations of a vortex 
patch. Dotted line: the Kirchoff ellipse. Full line: perturba- 
tive solution, (a) &2=0.2, (b) 62 = 0.3. 

In the case of deformations with higher angular de- 



pendences there are no analytic solutions beyond the lin- 
ear approximation. In this case, the angular frequen- 
cies are givena to zeroth-order in the deformation by 
f2 = (w— l)/2m, which coincides with the values for Qf 



[see Eqs 
Zabusk 



and 



B5 



Numerical solutions by Deem and 
for I = 3, 4 also compare favorably with the 
perturbative results, as can be seen in Fig. |^. In fact, 
the angular velocities determined from Table II coincide 
with those shown in Ref. |ll| to all significant figures of 
that paper. 



TABLE II. Stationary deformations of a vortex patch. Angular frequencies and lowest harmonics of the deformation for 
L — 2,3, 4, 5 [the equivalent results for the 2DES are summarized in Table |J. 



£7 = Q,/lu p 



64/ 



i - 61 + bi 



1 - 26| - 206| 



§ - 36! - 57bj 



4&! - 1526| 



3 h 2 1 h 4 

2°2 - 2°2 

§&! + 166| 
|6| + 1726| 



5 l3 
2°2 



■M . 
— 4 



231 i,4 



715 l4 
-9- "4 




FIG. 9. Vortex-patch "eigenstates" with L = 2> and L — 4. (a)-(d) Obtained numerically by Deem and Zabusky [figures are 
scanned from Ref. (e)-(h) Determined perturbatively. The lowest harmonic bh is taken in each case from Table I of Ref. 
[n]: 63 = 0.096, i>3 = 0.11, 64 = 0.048 and 64 = 0.070 respectively. The solutions are essentially identical and are too close to 
compare effectively in the same picture. 
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APPENDIX C: THE KIRCHHOFF ELLIPSE AS A 
SOLUTION OF A CONTOUR DYNAMICS 
PROBLEM 

It is illustrative to show, starting from the contour dy- 
namics, that an ellipse is indeed an invariant deformation 
for a vortex patch. Additionally it can be shown that this 
is an inherent property of the logarithmic kernel [Eq. B3] 
and that an ellipse it is not a solution for the 2DES. 

Consider the following parametrization of an ellipse: 

x — (a + b) cos 77 , 
y — (a — b) sin 77 . 

It is evident that the radius is given by 

(C2) 



(CI) 



r{rf) = \/ d 2 + 2ab cos 2ij + b 2 , 
and that tangential and normal vectors are given by 
dr 

e z x r (a — b) cos r/e x + (a + b) sin 77 e y 



- (a + b) sin -q e x + (a — b) cos 77 e y , (C3) 



where |t| = a 2 — 2ab cos 2i] + b 2 . The "rigid-body ro- 
tation" condition [Eq. (|l3[)1 can then be written as 



U = 



r(v)\ 



2ab sin 2rj . 



(C4) 



In addition the distance between two points on the ellipse 
takes the simple form 



J? = |r(T7) - r(^)| s 

-4(a 2 + fe 2 )sin 2 ^' "'' 



(C5) 



2ab 
i 2 + b 2 



cos(r/ + 77') 



and the dot product involved in Eqs. 
given by 



§ and (B4) is 



11(77) ' ^(77') 



1 



\r(v)\ 



{a 2 - b 2 ) sin(7/ - 77') . (C6) 



1. Vortex patches — exact solution 

It is now simple to show that an ellipse is, indeed, a uniformly rotating shape for the vortex patch [Eqs. 



(B4)], since the normal velocity is given by 



and 



Uy 



4vr|x(77)| J 



dr)' In 



4vr|x(77)| J 



2tt 

dx In 



4(a 2 + 6 2 )sin 2 (^-JL) 



2ab 



b 2 



cos (77 + rf) 



2ab 

1 5 cos x 

a 2 + b 2 



x (sin 2rj cos x — cos 2r\ sin x) 



sin (77 — 77') 



(C7) 
(C8) 



where in going from Eq. ( |C7| ) to (C8) we eliminated the 
term inside the first braces in the logarithm due to sym- 
metry and changed variables to x — 77 + 77'. Finally, the 
term proportional to cos 277 vanishes upon integration 
and we are left with a simple integral, proportional to 
sin 277. As long as a>b (for a — b the ellipse has collapsed 
into a line), the integral exists in closed form: 



(a 2 - b 2 )b 
2a\T{j])\ 



(C9) 



which, by direct comparison with Eq. (C4), yields the 
angular velocity 



v_ 4 I a 2 



(CIO) 



Consider now the parametrization given by Eq. (jlj) 
with bi given by the first row in Table II. The maximum 
and minimum radii correspond to R(ip) for ip equal to 
and 7r/2 respectively (for 62 > 0). It is easy to see that 
a~l+36|-6| and b^2b 2 +hb\, thus Q/o;p~l/2-&|+6|, 
as shown in Table II. A simple Fourier analysis of Eq. 
( |C2| ) also results in coefficients 6/ which agree with the 
perturbative solution. 



2. Two dimensional electron systems- 
solution 



-no exact 



It is also simple to see why an ellipse is not a station- 
ary solution of the two-dimensional electron system. In- 
stead of the logarithm, one has to deal with 1 / vR?, and 
the elimination of the first term inside the braces is not 
possible. The resulting integrands depend on 77 in a non- 
trivial way, and the normal velocity is not proportional 
to |t(?7)| -1 sin 277. 



APPENDIX D: GEOMETRY OF PLANAR 
CURVE MOTION 



For completeness, it is worth considering some general 
features of the planar curve motion.EZI Consider a curve 
described by some parametrization r(a), where a is a 
parameter defined on a fixed interval (see Fig. |l0|). It is 
then possible to consider tangent and normal unit vectors 
defined by 
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T<a) , , . dr(a) 

ts kWI» where T(a) = r " = ~toT ' (Dl) 

n = e r x i . 



The arc-length s(a) corresponds to the length of the 
curve from some arbitrary point to r(a) and is defined by 
ds = yj dx 2 + dy 2 = ^fgda, where the metric g is defined 
by g = t • t. We then have the Frenet-Serret relations 



r s = t , 
t.„ = -o, 
n, = « t , 



(D2) 



which define the curvature k. It is also possible to de- 
fine the curve in terms of its tangent angle 8(s), where 
8. = K. 




FIG. 10. Parametrization of a closed curve by a parame- 
ter a defined on a fixed interval. The illustration shows the 
arc-length s, the normal and tangential unit vectors n and t 
and the tangent angle 9. 

The kinematics of the curve can be determined if the 
velocity of the curve r t is determined. It is convenient 
to decompose this velocity into its normal and tangential 
velocities (at fixed a): 



r t = Un + Wt. 



(D3) 



In general, U and W may be arbitrary functionals of r(s) 
and its derivatives. It is easy to show that 



n t =-{U s -nW)t, 
tt = (U s - nW)n , 
(y/g)t = W a + y /gnU. 



(D4) 



The length C and area A of the curve are given by 

C= I ds = J^gda, (D5) 

A = ^ I(r X t) • ds = - J '(r x r Q ) • e z da , (D6) 

and it is evident from the expressions above that their 
time derivatives are equal to 



C t = 



A t 



<£(KU + W s )ds= I nUds, (D7) 



U(s)ds, 



(D8) 



where in the rightmost part of Eq. (D7) we assumed that 
W was periodic. 

Consider now the following identity 

d d did g t d , d d 



dt ds dt ^fg da 2g 3 / 2 da ds dt 



Using Eq. (D4), we obtain the commutator between the 
arc-length and time derivatives 



d_ d_ 

ds'di 



(W s + kU) ^ 



(D10) 



We are finally able to determine a kinematic equation 
for the curvature k. While the procedure is completely 
general, this becomes particularly important when the 
dynamics obeys a geometric law of motion, that is when 
U and W are functions of the curvature and its deriva- 



tives only. Using Eqs. (D4) and (D1C) to calculate the 
time derivative of the curvature at fixed a we find that 



dn 



dt 



= -[K 2 + d ss ]U + K S W . 



(Dll) 



It is convenient, however, to write a differential equa- 
tion in a parametrization-independent form, that is the 
time derivative should be evaluated at fixed arc-length s. 
Since <9 t | s = d t \ a — n s Jq[kU + W s >] ds' , we have 



dn 



dt 



-[k 2 + d ss ] U + k s W - k s ( [kU + W s >] ds' , 

Jo 



(D12) 



which is the form used in Eq. (| 



APPENDIX E: ELLIPTIC INTEGRALS, 
ELLIPTIC FUNCTIONS AND BOUNDARY 
CONDITIONS 

A quick inspection of Eq. fll^ ) reveals the following pos, 
sible alternative solutions, in terms of elliptic functions,E3 
which depend on the nature of the zeros of V(k): 

1. Four real zeros a>K>/3>"f>5 (Fig. ||-d): 



z(k) 



y/(a-i)(l3-6)' V("-/ 5 )( K -7) 
\l = (a-p){i-8)/{a- 1 ){f3-8), 

which can be inverted to read 

0(ja) + y(a - (3) sn 2 [c 4 z, A 4 



(7 — a) + (a — f3) sn 2 [04 z, A4 



(El) 
(E2) 

(E3) 



with cf = (a — 7)(/3 — S)/16. The period of k(z) 
is given in terms of the complete elliptic integral 
L 4 = (2/c 4 )K{\ 4 ). 
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2. Two real and two complex conjugate zeros a>re> 
(3, m±in (Fig. ||-a,b,c): 



z(k) = 



■F\2 tan 



m 

2 l \2 , 2 

p = (to — a) + n 

2 /„, o\2 . 2 



/ g(a— «) 



(to - py + n l 

>4 = [(Kmax - Kmin) 2 ~ (p ~ <?) 2 ] , 

which can be inverted to read 

(f3p + aq) + (/3p - aq) cn[y/pqz/2, A 2 ] 



(E4) 
(E5) 



k{z) = 



(p + q) + {p- q) cn[^/pqz/2, A 2 ] 



(E6) 



The period of re(z) is given in this case by by 
Li = (&/Jpq)K(\ 2 ). 



Since there is no cubic term in the potential V(k), the 
sum of all roots must vanish. That leaves only one free 
parameter in each case, once the minimum and maximum 
curvatures are fixed. As mentioned in Sec. VII , this free 



parameter is determined by the boundary condition that 
the curve is closed and without self-crossings. For the 
case of four real zeros, it is not possible to find solu- 
tions that satisfy these conditions. It is still possible to 
find beautiful closed curves (Fig. O), but these do not 
correspond to physical solutions for the problem under 
consideration. The case with two real and two complex 
conjugate solutions does have physical solutions and is 
discussed in detail in Sec. VII. 






FIG. 11. Examples of closed boundaries for the solutions 
corresponding to four real zeros of V(k). Note that all curves 
cross themselves at least once and are not physical solutions 
for the problem of interest. 



the curvature satisfied k(sA) — re(s — ct) were obtained 
by twice integrating Eq. (p^), thus resulting in Eq. Q), 
which can then be reduced to a simple problem in quadra- 
tures [Eq. (E^)] and the solutions for the curvature in 
terms of Jacobi elliptic functions were discussed in that 
context [see also Appendix || . 

In the case of an infinite curve, however, k and all its 
arc-length derivatives should vanish as s — > ±oo, and the 
constants of integration a = b = in Eqs. (thy) and (Eq), 



so that V(k) = k 4 /4- 
can then be found as 



ck . The traveling wave solutions 



k(z) = 2\[^c sech[\/— c is — ct)] , c < . 



(Fl) 



The tangent angle and the shape can then be obtained 
by a simple integration: 



z) = nds = 4arctan<tanh 



(s-ct) 



ly 



(s - ct) 



ds' 



■ tanh[^-2 £ (s - ct)] 



7T, (F2) 



(F3) 



Unfortunately these curves represent a small loop trav- 
eling along the boundary, and self-crossing solutions are 
not possible for the boundary of a physical system. 

It is, however, possible to have more complicated solu- 
tions that have a traveling envelope with time dependent 
oscillations within it. One such example is the "breather" 
solutionJi3 which loosely speaking, corresponds to a pair 
of bound solitons: 

, . d {I sm\k s - k(k 2 - 3l 2 )t] 

K (s,t) = -4- arctan j - -^—^-^ 

(F4) 

where I and k are arbitrary. The value for the tangent 
angle in this case is evident since 6 S = re; however, the 
shape of the curve requires numerical integration. Figure 
|l2] shows a sequence corresponding to the motion of one 
example of a breather. 




APPENDIX F: THE STRAIGHT INFINITE EDGE 



In Sec. VII we discussed the invariant shapes of a closed 
curve when the curvature satisfies a modified Korteweg- 
de Vries equation (E^). Single soliton solutions where 



FIG. 12. A breather soliton for k = 2 and 1 = 1. In this 
time sequence At = 0.125. 



14 



1 See for example H. Goldstein, Classical Mechanics 
(Addison-Wesley, Reading, 1980), 2nd. Ed., Chap. 1,2 and 
8. 

2 L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perga- 
mon Press, Oxford, 1959) §12. 

3 R. G. Holt and E. H. Trinh, Phys. Rev. Lett. 77, 1274 
(1996); R. E. Apfel et al, Phys. Rev. Lett. 78, 1912 (1997). 

4 A. Ludu and J. P. Draayer, Phys. Rev. Lett. 80, 2125 
(1998). 

5 T. M. O'Neil, Physics Today 52 vol. 2, 24 (1999), and ref- 
erences therein. 

6 A. Bohr, Mat. Fys. Medd. Dan. Vid. Selsk. 26 (1952) no. 
14; A. Bohr and B. R. Mottelson, Nuclear Structure, Vol. 
1 (Benjamin, New York, 1969). 

7 E. A. Overman and N. J. Zabusky, Phys. Rev. Lett. 45, 
1693 (1980). E.A. Overman, N.J. Zabusky and S.L. Os- 
sakow, Phys. Fluids 26, 1139 (1983). 

8 R. E. Rosensweig, Ferrohydrodynamics (Cambridge Univ. 
Press, 1985); S. A. Langer, R. E. Goldstein and D. P. Jack- 
son, Phys. Rev. A 46, 4894 (1992); A. T. Dorsey and R. 
E. Goldstein, Phys. Rev. B 57, 3058 (1998). 

9 Sir H. Lamb, Hydrodynamics (Dover, New York 1932) Sect. 
158 and 159. 

10 N. J. Zabusky, M. H. Hughes and K. V. Roberts, J. Comp. 
Phys. 30, 96 (1979); N. J. Zabusky and E. A. Overman, 
ibid 52, 351 (1983). 

11 G. S. Deem and N. J. Zabusky, Phys. Rev. Lett. 40, 859 
(1978). 

12 R. E. Goldstein and D. M. Petrich, Phys. Rev. Lett. 69, 
555 (1992). 

13 See X. G. Wen, Int. J. Mod. Phys. B 6, 1711 (1992); and C. 
L. Kane and M. P. A. Fisher, in Perspectives in Quantum 
Hall Effects: Novel Quantum Liquids in Low-Dimensional 
Semiconductor Structures, edited by S. Das Sarma and A. 
Pinczuk (Wiley, New York, 1997), pp. 109-159. 

14 C. Wexler and A. T. Dorsey, Phys. Rev. Lett. 82, 620 



(1999). 

15 D. J. Korteweg and G. de Vries, Phil. Mag. (5), 39, 422 
(1895); A. Das, Integrable Models (World Scientific, Sin- 
gapore, 1989); P. G. Drazin and R. S. Johnson, Solitons: 
an Introduction (Cambridge University Press, Cambridge, 
1989). 

16 R. E. Goldstein and D. M. Petrich, Phys. Rev. Lett. 67, 
3203 (1991); in Singularities in Fluids, Plasmas and Optics 
(Kluwer Acad., 1993), pp. 93-109. 

17 V. A. Volkov and S. A. Mikhailov, Zh. Eksp. Teor. Fiz. 94, 
217 (1988) [Sov. Phys. JETP 67, 1639 (1988)]. 

18 S. Giovanazzi, L. Pitaevskii, and S. Stringari, Phys. Rev. 
Lett. 72, 3230 (1994); A. Cappelli, C. A. Trugenberger, 
and G. R. Zemba, Ann. Phys. 246, 86 (1996). 

19 S. Iso, D. Karabali, and B. Sakita, Nucl. Phys. B 388, 700 
(1992); Phys. Lett. B 296, 143 (1992); D. Karabali, Nucl. 
Phys. B 428, 531 (1994). 

20 I. L. Aleiner and L. I. Glazman, Phys. Rev. Lett. 72, 2935 
(1994). I. L. Aleiner, D. Yue and L. I. Glazman, Phys. Rev. 
B 51, 13467 (1995). 

21 G. Ernst, R.J. Haug, J. Khul, K. von Klitzing and K. Eberl, 
Phys. Rev. Lett. 77, 4245 (1996). 

22 M. Abramowitz and I. Stegun, Handbook of Mathematical 
Functions, (Dover, New York, 1965). 

23 S. J. Allen, H. L. Stormer and J. C. M. Hwang, Phys. Rev. 
B 28, 4875 (1983); M. Wassermeier et al, Phys. Rev. B 
41, 10287 (1990); I. Grodnesky, D. Heitman and K. von 
Klitzing, Phys. Rev. Lett. 67, 1019 (1991). 

24 R. C. Ashoori H. L. Stormer, L. N. Pfeiffer, K. W. Baldwin, 
and K. West, Phys. Rev. B 45, 3894 (1992). 

25 G. Ernst, N. B. Zhitenev, R. J. Haug, and K. von Klitzing, 
Phys. Rev. Lett. 79, 3748 (1997). 

26 See G. K. Batchelor, An Introduction to Fluid Dynamics 
(Cambridge University Press, Cambridge, 1967), Ch. 7. 

27 R. C. Brower, D. A. Kessler, J. H. Koplik and J. H. Levine, 
Phys. Rev. A 29, 1335 (1984). 

28 J. Miller, Phys. Rev. Lett. 65, 2137 (1990); J. Miller, P. 
B. Weichman, and M. C. Cross, Phys. Rev. A 45, 2328 
(1992). 

29 N. B. Zhitenev, R. J. Haug, K. v. Klitzing, and K. Eberl, 
Phys. Rev. B 52, 11277 (1995). 



15 



